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In this work, we use a passive advection model for ion temperature fluctuations, in order to 
\Q , investigate their effects on Doppler Spectral line shapes. The relevance of the model is discussed 

in the framework of the Braginskii equations, and the subsequent Probability Density Function 
evaluation relies on results obtained in neutral fluids. The resulting Doppler line profiles are shown 
£SJ , to exhibit characteristic exponential tails. 
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I. INTRODUCTION 

The understanding of turbulence in fusion devices has made great progress since it was first realized that turbulent 



fluctuations could explain the anomalously high level of cross magnetic field transport which was plaguing experiments. 
Theoretical studies, many of which have focused on the identification of underlying linear or non linear instabilities pj, 
have produced valuable results. Direct Numerical Simulation, which is the only tool able to encompass the complexity 
of the problem (non linear equations, multiple fields and complex geometry) , has shed additional light on the respective 
contributions of these instabilities (e.g. 0)- Edge plasma turbulence has specific properties and is, in particular, 
O characterized by high fluctuation levels which can rise up to 30 % 0, 13 ■ The development of sophisticated diagnostics 
^ | methods now provides detailed experimental data on density fluctuations (e. g. [5J), which impose more stringent 
O ■ constraints on the theory. However,very few measurements of ion temperature fluctuations have been reported to 
date. A notable exception is Ref. 0, in which High Frequency Charge Exchange Recombination Spectroscopy is 
used to show that ion temperature fluctuations rates can be large. In previous works we have investigated the 

' possible role of turbulence on Doppler line shapes obtained by passive spectroscopy. Ion temperature fluctuations 
were shown to lead to significant modifications in the line wings, provided their PDF (Probability Density Function) 
had fat tails. This suggests that passive spectroscopy, for which the experimental set-up is relatively simple, might 
t-H ! provide information on ion temperature fluctuations. Here, making use of results obtained in neutral fluids, we show 
^ ■ that even a very simplified model for ion temperature fluctuations, i.e. passive advection, leads to non Gaussian 
PDFs. The corresponding asymptotic behavior of the Doppler line profiles is obtained analytically. 
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II. DOPPLER SPECTRA IN TURBULENT PLASMAS 

: 

The spectra / m (AA) measured using passive spectroscopy is obtained from integration both during the acquisition 
""^5 ' time of the spectrometer r and along the line of sight (LOS) (Oz). The corresponding expression is 

o : 

i f r i r 

I m (AX) = -j dt- J dzB(z,t)I(AX,z,t), (1) 

where L is the emitting zone length, B(z, t) the brightness of the line and /(AA, z, t) the local line shape. Both of these 
. . quantities are determined by the plasma parameters and the line under study. Here we consider D a (n = 3 to n = 2 
transition, Ao = 6561 A), which is intense and routinely monitored. The brightness essentially depends on the electron 
density and temperature. For ionizing edge plasma conditions, that is N e = 10 19 — 10 20 m~ 3 and T e = 10 — 100 eV, 
Stark effects can be safely neglected. We will therefore deal with the Doppler line profile of one Zeeman component. 
Here we only consider the contribution from the class of neutrals created locally by charge exchange, I(AX, z,t) — > 
J CX (AA,X). In fact, neutrals created via other channels,i.e molecular dissociation and charge exchange in the inner 
plasma are not strongly coupled to the ions at the edge jg, i.e. to edge turbulence. These neutrals can be accounted 
for using neutral transport codes. The charge exchange neutrals Velocity Distribution Function (VDF) along the 
LOS, f cx (v,X), is related to the Doppler profile by f cx (v,X.)dv — I cx (AX,X)dAX. Here it has been made explicit 
that the space and time dependence of the spectra and VDF occurs only through the local ion fluid fields X(,z,t) 
(local number and charge densities, species temperatures, and flow velocities), since the ion VDF is assumed to be 
a local Maxwellian. We consider the case where the LOS is such that the average gradients of these fields along its 



2 



direction (e.g., parallel to the magnetic field lines) are weak . It was shown in Ref. [aQl that the measured spectra 
can then be written as 



I™(A\) = / dXW{X)B{X)r x {A\,X), (2) 

where VF(X) is the joint PDF of the fields. To proceed, the role of the different fluid fields can be investigated 
separately. For Doppler line shapes, fluid velocity and temperature fluctuations give rise to the main effects, because 
the shape of the local VDF f cx does not depend on N e . We further simplify the problem by considering ion temperature 
fluctuations alone. The rationale behind that assumption is twofold : first, the role of velocity fluctuations is well 
known in plasma spectroscopy . Secondly, the analysis carried out in our previous work showed that line wings are 
significantly affected by temperature fluctuations characterized by a slowly decaying tail. Therefore, we expect that 
inclusion of velocity fluctuations would not significatively shroud the conclusions drawn here, since line wings could 
only be affected by fluid velocities several times larger than the ion thermal velocity. Moreover, if the LOS is parallel 
to the magnetic fields, there are conditions under which parallel velocity fluctuations are small while temperature 
fluctuations are not As a result, the remainder of this contribution will be devoted to temperature fluctuations. 

III. SIMPLIFIED MODEL FOR THE ION TEMPERATURE 

We start from the Braginskii equations an d focus on the ion thermal balance 

^dtimTi) + -n;T;V • v, + 7T :: a + V • q, = Q l {T l ,v), (3) 

which we simplify using the drift wave ordering. Here d t = d t + v; • V stands for the convective derivative, q^ for the 
thermal flux, 7r Q ^ is the stress tensor, a a p = d a vp and Qi is the energy exchange with other species. We consider 
fluctuations having a small scale / with respect to the quasi equilibrium background of radial gradient length Lj_. 
Typically p s < I, where p s is the ion Larmor radius calculated with the electron temperature. For each field X, we 
distinguish between its time averaged part X and its fluctuating part X. The small scale parameter with respect to 
which the fluid equations are expended is 8 = l/L±. The perpendicular velocity is given by Vi± — ve + v* + vf , that 
is by the sum of the electric drift, ion diamagnetic and ion polarization velocities. The diamagnetic velocity does not 
directly advect temperature, because for low (3 (the ratio of the kinetic to the magnetic pressure) its effect cancels out 
with that of the diamagnetic thermal flux, up to magnetic curvature terms which will be neglected in the following. 
The same approximations ensure that V ■ ~ 0, i.e. the fluctuating electric drift velocity flow is incompressible. 
The ion polarization velocity is of higher order in 8 and we further neglect parallel advection, which is justified in the 
drift ordering |fl| . The resulting equation for ion temperature fluctuations is then 

[d t + (v E + v E ) ■ V - k 1± V 2 ] Ti = S(r, t), (4) 

where Ki± is the ion perpendicular thermal diffusivity, and all the terms on the r.h.s. have been lumped into a source, 
that is S^r, t) — ve ■ VT; — ve • VTi + Qi. The energy exchanges between electrons and ions Qi e and ions and neutrals 
Qi n are neglected on the turbulent time scales, i.e. Qi = Qi — Qi — 0. Also, the second term is zero by homogeneity, so 
that S(r, t) =Ve' VT$. The mean temperature gradient therefore naturally provides a source term for the fluctuations 
through its coupling to v^. Since the scale of the turbulent fluctuations is such that I <C L±, in the following the 
gradient will be assumed to be constant and directed along x, namely V x Ti — g. Advection by the mean flow ve in 
Eq. can be formally removed by introducing suitable Lagrangian coordinates. This does not interfere with the 
following developments, so that we will assume v^; = 0. Rigorously speaking, the evolution of v^; is not independent 
from that of Ti (see for example the system of equations considered in Ref. 0). However, the ion temperature does 
not play any role in the parallel electron dynamics, in contrast to the electron temperature which is strongly coupled 
to the electrostatic potential, and hence to v^. In addition, many if not most edge turbulence studies carried out so 
far did not include ion temperature fluctuations, by considering cold ions (Ti <c T e ). This approximation is justified 
partly by the fact that ion temperature does not play a fundamental role for the instabilities which are thought to be 
dominant in edge plasmas (collisional drift waves, resistive ballooning). In fact, inclusion of Ti does not qualitatively 
change the nature of the turbulence observed in the simulations 0, and the few experimental results available are 
not inconsistent with passive advection, as pointed out in 4j. In this work, we therefore treat ion temperature as 
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a passive scalar driven by the fluctuating electric drift velocity field. This has the major advantage to lead to an 
analytically tractable model. Let us now specify the typical values of the Peclet and the Prandtl numbers, which 
will play a role in the next section, for edge plasma turbulence. The Peclet number, Pe = lv/Ki±, where I and v are 
the typical length scales and velocity of the turbulent fluctuations, controls the relative importance of advection and 
diffusion. If we take I ~ 10p s ~ 2 x I0~ 3 m and v ~ ve ~ v* ~ Ti/(eBL±) ~ 75 m.s -1 with L± = 2 cm, we get 
Pe ~ 10 2 3> 1, so that advection dominates the heat transport. The Prandtl number Pr = Vi±/Ki±, where Ui± is the 
perpendicular ion viscosity is of order unity since Vi±/ni± ~ 1 [T(i| . 

IV. REVIEW OF ANALYTICAL CALCULATIONS OF THE PASSIVE SCALAR PDF 

Since and v are rapidly fluctuating quantities, we adopt a stochastic description of the velocity field. For the 
sake of tractability, the latter is described by a Gaussian probability density functional, whose correlation function is 
chosen so as to reproduce some of the main features of turbulence [I3 , that is 

C a(3 (r,t) = (v a {T,t)vp(fl,Q)) = n a/3 7£exp(-|r|/£)tf(i) (5) 

where U a p(r) ensures consistency with incompressibihty [lflj (at zeroth order in r, H a b = S a b), and the white noise 
limit has been taken. Here, the correlation length £ and the typical velocity V are such that £ ~ I and V ~ v. The 
time average is thus replaced by an ensemble average denoted by (•). In this section, we show how the ion temperature 
PDF W(9) = (6(9 — Tj(r, <))), or equivalently its Fourier transform Z{\) = (exp(— iATj(r, £))), can be obtained from 
Eq. (@J. Here, 9 is a sample space variable, and A its conjugated Fourier variable. Mathematically, the brackets 
denote functional integration over the field v at all space and time points. The PDF does not depend on t and r 
because of stationarity and homogeneity, the latter stemming from the fact that the average temperature gradient is 
constant. The derivation presented here is drawn from the neutral fluid community, which has given much attention 
to the passive advection problem. Indeed, its study allows to unravel important effects also at play in Navier-Stokes 
turbulence. The passive scalar equation is linear in X^, but bilinear in the fluctuations and can thus generate non 
linear effects. The scalar PDF W(9) can therefore have strong non Gaussian features even though the velocity field is 
Gaussian. This is very well illustrated by the early work of Sinai and Yakhot, who found power law tails for the PDF 
in homogeneous decaying turbulence |l4j| . The derivation relies on a transport equation for the PDF, which involves 
coefficients whose temperature dependence has to be determined from rather uncontrolled assumptions. This difficulty 
arises because diffusion introduces correlations between the scalar field and its gradients. The physical reason for this 
is most easily understood using a Lagrangian picture. In the weak thermal diffusion limit (i.e. large Peclet number 
Pe 3> 1, see below), advection is the dominant initial effect and the fluid particles follow the Lagrangian trajectories 
defined by d t R(t) = v(R(t),f|Ro) where R(t = 0) = Ro = r. The velocity field is linearized around the Lagrangian 
trajectory (v a (r,t) — > v a (R(t), t) + o- a p(TL(t), t) (rp — Rp(t))), where <j a/ g(R(f), t) — d a vp(R(t),t) is the strain field 
along the trajectory. Two points in a fluid element, initially close will be displaced by this strain field. Due to the 
condition of incompressibihty, one direction will involve an exponentially large separation as a function of time. As a 
consequence of this so-called " advective stretching" , the spectra of temperature fluctuations are extended to smaller 
and smaller scales. When the diffusive scale is eventually reached, regions of large gradient are selectively dissipated, 
introducing correlations between the temperature field and its gradient [15| . The PDF obtained in Ref. pjj has been 
used as an example in our previous work 0, 01 1 but from the section 3 it is clear that an additive noise should also be 
taken into account. The corresponding transport equation for the PDF could straightforwardly be written down, but 
finding its solution would require further assumptions on the correlations between Ti and S [Lsl ] . From the discussion 
above it is obvious that a Lagrangian description would allow to capture more of the physics. The remainder of this 
section will be devoted to show how this unfolds. The first difficulty when carrying out the average in Z(\) lies in the 
fact that there are two different sources of noise associated with ve, one additive (S — gvEx) and one multiplicative 
(ve ■ V). An elegant way to disentangle their respective contributions is to calculate the Green's function G for Eq. 
0, satisfying 

Ti(i,t)= f dr'dt'G(r,t\r',t')gvE X (r',t'). (6) 

This has the advantage to allow averaging on the additive noise as a first step. The latter is independent of the 
statistics of ve in G(r, t\r', t') in the limit of white noise. Physically speaking this occurs when the time scale for 
velocity field correlations is short compared to other relevant time scales for Tj. In the following, all results quoted 



4 



are in this limit. The resulting expression for the characteristic function is 



Z{\) = (^cxp \-X 2 J dt J dkD(k)|G(k,t)rjy , (7) 

where {-)m stand for the ensemble average over the multiplicative noise, i.e. the v dependance of the Green function. 
D(k) is the Fourier transform of the additive noise correlation function D(r) — (S(0,i)S(r,t))A = g 2 C xx (r). This 
approach was first carried out in Ref. ^3 using a path integral formulation to calculate the Green's function, assuming 
the above linearization of the velocity field around the Lagrangian tr ajec tory. Strictly speaking, the latter is valid for 
Pr 3> 1, but is thought to be also relevant to the cases where Pr ~ 1 |X3| - The same approximation allows to expand 
C xx {v) in terms of |r|/£. The Green's function, Fourier transformed with respect to the starting point, is then found 
to be 

G(0, 0|k , t) = e *a-fc°.°) exp (-^L J* dt'k 2 (t')\ , (8) 
where the wave number k(i) is defined by the time ordered product (7~ stands for latest times on left) 



k a (t) = T 



exp (- £ a(t')dt' 



hop. (9) 

a/3 



Here r(t; 0, 0) is the backwards Lagrangian trajectory, i.e. the position at which the fluid particle which is at time t = 
at position r = was at time t < 0, not taking into account diffusion. The second term in Eq. J3J) is a damping term 
originating from diffusion, where the time dependence of k(i) describes the effect of stretching integrated along the 
Lagrangian trajectory. The norm of k grows exponentially for the velocity field considered here, so that this second 
term basically behaves as a step function at a time t* , which can be defined as |G(0, 0|ko, t*)\ ~ 1/2. The larger the 
strain along the trajectory, the smaller t*. For t < t* the fluid particle follow the diffusionless Lagrangian trajectory 
with constant temperature, and then the latter is quickly homogenized because of the large gradients induced by 
stretching. Since the modulus of the Green's function only depends on the velocity field through the strain tensor 
a (see Eq. the velocity average is then converted into an average over realizations of a. Carrying out the 

calculation is tedious, but only involves some additional mild approximations [l2f . The resulting PDF has a Gaussian 
core and an asymptotical behavior which is exponential, that is using non dimensional quantities 

WW-^e-rW (10) 

where 7 = (d + 6)/2d is the secular growth rate, and d the dimensionality of the problem (7 = 2 in 2D) The 
same result has been obtained using Jensen's field theoretical formalism [Vq. Taking an additive noise into account, 
as required by the derivation of Eq. (@J , therefore leads to a PDF which has an exponential tail instead of the power 
law tail obtained for decaying turbulence. The PDF is however still markedly non Gaussian. 

V. CALCULATIONS OF DOPPLER SPECTRA 

We now consider the modifications resulting from temperature fluctuations on line shapes, given the PDF obtained 
in the previous section. The Doppler line shape / m (AA) can be calculated numerically from Eq. © for any given 
W(9), noting that the brightness of the line B does not sensitively depend on the ion temperature. However, it is 
interesting first to obtain an analytical expression of the profile. This can be done for PDFs behaving asymptotically 
as W{6) ~ 6~ a exp — (9/9q)^ . The asymptotic behavior of the measured Doppler profile can then be calculated using 
a saddle point approximation, which leads to 

7 m (AA)cx l__exp(-C(0 o ,/3)|AA|^T) , (11) 

|AA|^r^ V ) 
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where C{6q,0) = ((3 + l)/(/3C0o)^ /3+1 , with £ = A 2 ,/ '(c 2 mi), rrii being the ion mass and c the speed of light. For 
(3 = 0, the measured profile decays algebraically in accordance with Ref. 0,0) while for an exponential PDF ((3 = 1), 
its decay is exponential too. This implies that ion temperature fluctuations characterized by an exponential PDF 
would lead to a conspicuous behavior in the line wings. As an example, we consider the following temperature PDF 
(adapted from Eq. (4.10) in Ref. [13), having a Gaussian core and exponential tails 



N I 9 f) 2 \ 
W(0) = ===== -(0 + T ) exp , (12) 

(gd)V9t(lC/VO^Pe+\0\ 3/2 V 9Zg&>C/VOhiPc+y/Tcm\0\J 

where TV is a normalization constant, £ and V are defined by Eq. (JjJJ, and 1C = / +O ° C xx (R.(t), t)dt is the eddy 
diffusivity, where R(0) = 0. The effective Peclet number V£//C is such that JC/V£ = 1 according to Eq. JSJ. This 
VDF is plotted on Fig. for T = 50 eV, Pe = 10 2 , and two gradient values g£ = 10 eV and 15 eV. The resulting 
Dopplcr profiles are plotted on Fig. together with the Gaussian at 50 eV corresponding to the fluctuation free 
case. The center of the line is not strongly affected, whereas significant deviations occur in the line wings. The 
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FIG. 1: a) Plot of the PDF of Eq. JTJ for T = 50 eV, g£ = 10 eV, Pe=10 2 , V£/tC = 1 (dotted line), and g£ = 15 eV 
(dashed line) . b) Plot of the resulting Doppler profile on a logarithmic scale (solid and dashed lines), showing the asymptotic 
exponential behavior. The solid line corresponds to the Gaussian Doppler profile obtained for 50 eV. 



linear dependance characterizing an exponential fall off in log-linear scale is clearly seen. These results show that 
Doppler line profiles would indeed be affected by ion temperature fluctuations stemming from passive advection. 
Comparison to experimental spectra requires specifically designed measurements with a large dynamic range. With 
existing techniques, spectra can be recorded over five orders of magnitude using CCD detectors operated in specifically 
optimized readout mode |l8j| . In Ref. |l9j |. we presented spectra measured in the Tore Supra tokamak, the far line 
wings of which were behaving in a way consistent with a power law. However, reliable determination of the exact 
nature of the these deviations (exponential or power law behavior) would require improved measurements. In fact, 
the latter spectra were obtained as a by-product of routine measurements, and therefore line wings were recorded only 
over a limited range. These spectra nevertheless provided clear indications of the presence of an unexpected behavior 
in the far line wings, and are a strong incentive for further investigations. 



VI. CONCLUSIONS AND PERSPECTIVES 



We have shown that as a first approximation the ion temperature field in edge plasmas can be described as a 
passive scalar, advected by the fluctuating electric drift velocity field. The average temperature gradient naturally 
provides a forcing term for the fluctuations. The PDF of the passive scalar can be calculated using results obtained 
in neutral fluids, and have an exponential fall off. The latter translates into an exponential fall off for Doppler line 
wings. Therefore, passive spectroscopy might provide information on the statistical properties of ion temperature 
fluctuations, about which very few is known at this time. The effects on Doppler line shapes are observable when the 
temperature fluctuation PDF deviates significantly from gaussianity, so that their study could either uncover, or rule 
out such behavior. This work also provides a strong incentive to study the statistical properties of ion temperature 
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fluctuations by direct numerical simulation. This would allow to compare the PDF obtained numerically to those 
provided by our simplified model. 
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